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We calculate the antiproton flux due to relic neutralino annihilations, in a two-dimensional diffu- 
sion model compatible with stable and radioactive cosmic ray nuclei. We find that the uncertainty 
in the primary flux induced by the propagation parameters alone is about two orders of magnitude 
at low energies, and it is mainly determined by the lack of knowledge on the thickness of the diffusive 
halo. On the contrary, different dark matter density profiles do not significantly alter the flux: a 
NFW distribution produces fluxes which are at most 20% higher than an isothermal sphere. The 
most conservative choice for propagation parameters and dark matter distribution normalization, 
together with current data on cosmic antiprotons, cannot lead to any definitive constraint on the 
supersymmetric parameter space, neither in a low-energy effective MSSM, or in a minimal SUGRA 
scheme. However, if the best choice for propagation parameters - corresponding to a diffusive halo 
of L — 4 kpc - is adopted, some supersymmetric configurations with the neutralino mass m x < 100 
GeV should be considered as excluded. An enhancement flux factor - due for instance to a clumpy 
dark halo or to a higher local dark matter density - would imply a more severe cut on the super- 
symmetric parameters. 

PACS numbers: 95.35.+d,98.35.Gi,98.35.Pr,96.40.-z,98.70.Sa,11.30.Pb,12.60.Jv,95.30.Cq 



I. INTRODUCTION 

The recent WMAP measurements of the Cosmic Mi- 
crowave Background (CMB) anisotropics pj point to- 
wards a flat universe with a fraction J7a — 0.7 of the clo- 
sure density in the form of a negative pressure component 
- such as a cosmological constant or a scalar field - while 
the remaining Q m ~ 0.3 is matter. These conclusions are 
independently reached from the determination of the re- 
lation between the luminosity distance and the rcdshift 
of supernovae SNela 0] on the one hand and from the 
large scale structure (LSS) information from Galaxy and 
cluster surveys 0- The WMAP values of f2 m = 0.27 ±0.04 
and S!b = 0.044 ± 0.004 indicate that most of the matter 
is non-baryonic. The amount of baryonic matter J7b de- 
duced from CMB is in perfect agreement with the results 
from primordial nucleosynthesis and observations of the 
deuterium abundance in quasar absorption lines 

The nature of this astronomical dark matter has been 
challenging physicists for several decades and is still un- 
resolved. The favoured candidate is a weakly-interacting 
massive particle (WIMP). The so-called neutralino natu- 
rally arises in the framework of supersymmetric theories 
as the lightest combination of neutral higgsinos and gaug- 
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inos. Large efforts have been devoted to pin down these 
evading species @. Experimental techniques 0, 0, 111 
have been devised in order to be sensitive to the recoil 
energy which a neutralino may deposit as it crosses a ter- 
restrial detector. The annihilation photons from the neu- 
tralinos that populate the Milky- Way halo or extra- 
Galactic systems 01 are under scrutiny. As a matter 
of fact, a gamma-ray excess has been recently reported 
by HEGRA in the direction of the giant elliptical M87 
Antimatter cosmic-ray particles are also expected 
from neutralino annihilations inside our Galaxy. A subtle 
feature in the positron spectrum has actually been mea- 
sured by the heat Collaboration ^| f° r energies beyond 
7 GeV. 

This work is devoted to cosmic-ray antiprotons whose 
energy spectrum has already been measured with some 
accuracy. Much larger statistics will soon be collected by 
the AMS collaboration on board the ISS by the BESS-Polar 
long duration, balloon experiment and by the PAMELA 
satellite. Secondary antiprotons are naturally produced 
by the spallation of primary nuclei - mostly cosmic-ray 
protons and helions - on the diffuse gas of the Milky- 
Way ridge. If ncutralinos pervade our Galaxy, a pri- 
mary component adds up to that secondary distribution. 
The spectral distortion that ensues is expected a priori in 
the low-energy region for mere kinematic reasons fl3| : 
unlike for a neutralino annihilation, the center-of-mass 
frame of a spallation event is not at rest with respect to 
the Galaxy. In principle, an excess of low-energy antipro- 
tons is the signature of an unconventional production - 
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cither neutralino annihilation or small black holes evap- 
oration 01 f° r instance. However, because antiprotons 
undergo inelastic yet non-annihilating collisions with the 
interstellar material, the high-energy particles tend to 
lose energy and to populate the low-ener gy t ail of the 
spectrum that consequently is much flatter [15| than pre- 
viously estimated. This motivated the search of other 
cosmic-ray signatures such as antidcuterons 

[IE 03 An- 

tiproton production from primary cosmic-ray spallations 
is the natural background to any unconventional excess 
that would signal for instance the presence of the pu- 
tative neutralinos. The detailed calculation of that sec- 
ondary component [l8[ has required the determination 
of the propagation-diffusion parameters that are con- 
sistent with the B/C data By varying those pa- 
rameters over the entire range allowed by the cosmic-ray 
nuclei measurements, the theoretical uncertainty on the 
antiproton secondary flux has been found to be 9% from 
100 MeV to 1 GeV. It reaches a maximum of 24% at 10 
GeV and decreases to 10% at 100 GeV. This small scatter 
in the secondary antiproton spectrum is not surprising. 
Cosmic-ray nuclei such as LiBeB and secondary antipro- 
tons are both manufactured in the same place - the in- 
terstellar gas of the Galactic disk - through the same 
production mechanism - the spallation of primaries. 

The aim of this article is to calculate the supersym- 
metric cosmic-ray antiproton flux that arise from the 
diffusion-propagation parameter space and to estimate 
the uncertainties due to its spread. Since neutralinos an- 
nihilate all over the Milky Way and are not confined to 
the disk alone, we anticipate that the uncertainty in that 
primary component will be much larger than for secon- 
daries. 

The discusion will be shared in two main directions, 
brought to the fore by the structure of the equation de- 
scribing the primary flux. Production and propagation 
may be disentangled in the limit where the energy does 
not change as antiprotons travel. That is not strictly cor- 
rect as diffusive reacceleration as well as adiabatic and 
Coulomb losses generate a diffusion in energy space that 
is dicussed in Sec. IB II The elementary process of su- 
persymmetric antiproton production through neutralino 
annihilation is discussed in Sec. [H] both in an effective 
MSSM and in a supergravity-inspircd model. The two- 
zone propagation-diffusion model and the dependence of 
the primary antiproton flux $™ by on the propagation pa- 
rameters is described in Sec. IIIII The thickness L of 
the magnetic halo is naively expected to be the dom- 
inant source of uncertainty for <j>^ usy as the larger the 
confinement layers, the larger the fiducial volume where 
neutralino annihilations take place and the larger the su- 
persymmetric antiproton flux. Actually, L is combined 
with the diffusion coefficient K(E) and the Galactic wind 
velocity V c in order to get a precise value for the B/C 
ratio - and for the antiproton flux. We present the re- 
suls for the primary flux in Sec. IIVI where we estimate 
the uncertainties induced by the spread of the diffusion- 
propagation parameter space. We also briefly discuss the 



modifications of &p usy due to different choices in the dark 
matter distribution function, in its normalization and in 
the core radius values. In Sec. E] the comparison of 
the latest antiproton measurements with the antiproton 
fluxes predicted in different supersymmetric schemes will 
be discussed as a function of the propagation-diffusion 
parameters and of the neutralino Galactic distribution. 
Conclusions and perspectives will be presented in Sec. lVIl 

II. THE NEUTRALINO— INDUCED 
ANTIPROTONS: THE SOURCE TERM 

Antiprotons can be produced by self-annihilation of 
neutralinos in the Galactic halo. Dark matter neutralinos 
may be considered almost at rest in the Galactic frame 
since their average velocity is of the order of 300 km s" 1 . 
They are therefore highly non-relativistic. The produc- 
tion differential rate per unit volume and time is a func- 
tion of space coordinates (r,z defined in the Galactic rest 
frame) and antiproton kinetic energy Tp. It is defined as: 

^{r,z,T 9 ) = <w,) fl(T p ) (^^) 2 . (!) 

where {cr ann v)o denotes the average over the Galactic ve- 
locity distribution function of the neutralino pair annihi- 
lation cross section cr ann multiplied by the relative veloc- 
ity v, m x is the neutralino mass and p x {r, z) is the mass 
distribution function of neutralinos inside the Galactic 
halo. Since relic neutralinos behave as cold dark matter, 
their distribution has to follow the matter density profile 
Pdm(,t, z) of the Galactic halo: 

Px( r > z ) = £/?DM(r, z) (2) 

where £ parameterizes the fact that the dark halo may 
not be totally made of relic neutralinos (£ < 1). This 
would be the case when neutralinos are not responsible 
for the total amount of dark matter in the Universe, i.e. 
when their relic abundance Q, x h? is much smaller than 
the measured value for f^DM^ 2 - This is a situation which 
occurs in many supersymmetric models. It is reasonable 
to assume that £ has no space dependence and that is 
related to the relative amount of fl x h 2 with respect to 
^dm^ 2 - We will assume the standard definition: 

£ = min(l, CL x h 2 /0M) (3) 

where we have considered that neutralinos with £l x h < 
0.05 cannot be the dominant dark matter component. 

Finally, the second term in Eq. QJ, g(Tp) 7 denotes the 
antiproton differential spectrum per annihilation event, 
defined as: 

g( T p) = -w, 

0"ann dip 

= £BR( X ^F)(^j, (4) 
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FIG. 1: Scatter plot of the supersymmetric flux factor T = £ 2 (a" an nv}o/m-x as a function of the neutralino mass m x , calculated 
in the eMSSM. Panel (a) refers to supersymmetric configurations with the neutralino as a dominant dark matter component (i.e. 
0.05 < Q x h 2 < 0.3, and therefore a rescaling factor £ = 1). The light (green) circles show the eMSSM configurations for which 
the neutralino relic abundance lies in the preferred range for CDM, as determined by the combined WMAP+2dFGRS+Lyman-a 
analysis: 0.095 < flcDAih 2 < 0.131 Panel (b) refers to the neutralino as a subdominant dark matter particle (Q x h 2 < 0.05). 



where F lists the XX annihilation final-state particles 
which can subsequently produce antiprotons either di- 
rectly (hadronization when F = quarks or gluons, or 
through subsequent decay of F into quarks or gluons), 
BR(xX F) is the branching ratio for the production 
of F and dNp / dTp denotes the differential energy distri- 
bution of the antiprotons generated by F . For details on 
the calculation of g(Tp), see App. lAland Ref. pp| . 

The source term qp usy (r, z,Tp) is therefore a combina- 
tion of astrophysical factors (the dark matter density pro- 
file of the Galactic halo) and of particle physics proper- 
tics (the neutralino self annihilation cross section and the 
hadronization into antiprotons of the neutralino annihi- 
lation products). The astrophysical and particle physics 
quantities are factored out and can be studied separately. 
With the definitions given above, we can rewrite the an- 
tiproton source term as: 

q S p USy (r, z, T p ) = Tg(T p )pl M (r, z) (5) 

where we have defined the supersymmetric flux factor T 
as: 

t = e {<Tan, f° (6) 

which entirely depends on properties of supersymmetric 
models. 

We move now to discuss each term separately. 



A. The Galactic distribution of dark matter 

For most of our discussion, we will assume that the 
dark matter density distribution is described by a cored 
isothermal sphere. In terms of the radial distance r in 
the Galactic plane and of the vertical coordinate z, the 
density profile is: 

, s a 2 + R 2 Q 

PDM(r, z) = pi — (7) 

a + r + z z 

where a denotes the core radius of the dark halo and R Q 
is the distance of the Sun from the Galactic center. We 
have set a= 3.5 kpc and the IAU-recommcnded value 
Rq = 8.5 kpc. The value pi for the total local dark 
matter density is determined by taking into account the 
contribution given by the matter density of Eq. J7J to 
the local rotational velocity v TO t H3- The value of pi 
compatible with observations ranges from 0.18 GeV cm 
(for a low value of the rotational velocity, w ro t = 170 km 
s _1 and a non-maximal dark halo) to 0.71 GeV cm -3 (for 
V ro t = 270 km s _1 and a maximal dark halo) |2l|. The 
interval relative to the preferred value for the rotational 
velocity (v Tot = 220 km s" 1 ) is 0.30 GeV cm~ 3 < pi < 
0.47 GeV cm~ 3 j^. Our results will be presented for 
pi = 0.3 GeV cm" 3 . Since in the primary antiproton 
flux p 2 enters as a pure normalization factor, the fluxes 
obtained for different values of pi are easily rescaled. For 
instance, for pi = 0.47 GeV cm -3 , the antiproton fluxes 
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FIG. 2: The same as in Fig. calculated in the mSUGRA scheme. Panels (a) and (c) refer to cosmologically dominant 
neutralinos (0.05 < Q x h 2 < 0.3); panels (b) and (d) to subdominant neutralinos (Q x h 2 < 0.05). The upper row (panels (a) 
and (b)) is obtained for the universal soft -scalar mass mo smaller than 1 TeV (for these models, the neutralino is mostly a 
bino state); the lower row (panels (c) and (d)) refers to values of mo in excess of 1 TeV (in this case the neutralino may have 
a substantial higgsino component). The light (green) circles in panel (a) and (c) show the mSUGRA configurations for which 
the neutralino relic abundance lies in the preferred range for CDM, as determined by the combined WMAP+2dFGRS+Lyman-a 
analysis: 0.095 < ttcDAih 2 < 0.131 Q. 



would be a factor 2.45 higher than the corresponding ones B. Supersymmetric models 

for pi = 0.3 GeV cm- 3 . 

The existence of a relic particle in supersymmetric the- 
ories arises from the conservation of a symmetry, R- 
parity, which prevents the lightest of all the superpart- 
ners from decaying. The nature and the properties of 
We will come back to the topic of the dark matter this particle depend on the way supersymmetry is bro- 
density profile at the end of the paper, in Sec. IIV Bl ken. The neutralino can be the dark matter candidate in 
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FIG. 3: Antiproton differential energy distribution for pure annihilation channels as a function of the reduced kinetic energy 
xp = Tp/m x . Panel (a) refers to annihilation into a bb pair, for neutralino masses of: m x = 10,60, 100,300,500, 1000 GeV 
(from bottom to top); panel (b) refers to annihilation into a ZZ pair, for m x — 100, 300, 500, 1000 GeV (from top to bottom); 
panel (c) refers to annihilation into a scalar+pseudoscalar higgs pair hA, for m x = 300, 500, 1000 GeV (from top to bottom), 
and for: win = 120 GeV, rriA = 200 GeV, tan/? = 10 (ratio of higgs vev's) and a = (higgs mixing parameter); panel (d) refers 
to annihilation into a hZ pair, for m x — 300, 500, 1000 GeV (from top to bottom), and for: raj, = 120 GeV, tan/9 = 10 and 
a = 0. 



models where supersymmetry is broken through gravity- 
(or anomaly-) mediated mechanisms. The actual imple- 
mentation of a specific supersymmetric scheme depends 
on a number of assumptions on the structure of the model 
and on the relations among its parameters. This induces 
a large variability of the phenomenology of neutralino 
dark matter. In this paper we will consider neutralino 
dark matter in two different supersymmetric schemes: a 



low-energy effective-theory implementation of the Min- 
imal Supersymmetric Standard Model (eMSSM) and a 
minimal Supergravity model (mSUGRA) . 

The eMSSM is defined as an implementation of su- 
persymmetry directly at the electroweak scale, which is 
where the phenomenology of neutralino dark matter is 
actually studied. The large number of free parameters 
is reduced by a set of assumptions which are sufficient 
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FIG. 4: Branching ratios for the neutralino self-annihilation cross section in the eMSSM. Panel (a) shows the amount of the 
branching ratio for the annihilation into a fermion-antifermion final state (xx ~ * //)• Panels (b), (c) and (d) show the amount, 
relative to the // final state, of the annihilation into higgs bosons, gauge bosons and the mixed higgs-gauge bosons final state. 
Dark (red) points denote configuration with 0.05 < Q x h 2 < 0.3 (dominant relic neutralinos) . Light (green) circles indicate 
configuration with Q x h 2 < 0.05 (subdominant relic neutralinos). 



to shape the properties of the model at the electroweak 
scale. All the relevant parameters, which set the mass 
scales and couplings of all the supersymmetric particles 
(and of the higgs sector) are taken into account. The 
free parameters are: the gaugino mass parameter M 2 , 
the higgs mixing parameters /i, the ratio of the two higgs 
vacuum expectation values tan/3, the mass of the pseu- 
doscalar higgs rriA, a common soft -scalar mass for the 
squarks rriq, a common soft scalar mass for the sleptons 
to ; ~, a common dimcnsionlcss trilinear-parameter A for 



the third family (Ar — Aj = Arriq and Af = Amf, the 
trilincar parameters for the other families arc set equal 
to zero). We assume the standard grand unification re- 
lation between the [7(1) and the SU{2) gaugino mass 
parameters: Mi — 5/3 tan 2 Q-^Mi- The parameters will 
be varied in the following intervals: 100 GeV < M 2 < 
1000 GeV, 100 GeV < < 1000 GeV, 100 GeV < 
m A < 1000 GeV, 100 GeV < m^mj < 3000 GeV, 
1 < tan (3 < 50 and -3 < A < 3. 

A different approach is to embed supcrsymmctry in a 
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FIG. 5: The same as in Fig. |1J calculated in the mSUGRA scheme. Dark (red) points denote configuration with 0.05 < Q x h 2 < 
0.3 (dominant relic neutralinos) . Light (green) circles indicate configuration with fl x h < 0.05 (subdominant relic neutralinos) . 
Crosses (in blue) indicate the mSUGRA configurations with mo > 1 TeV. 



supergravity scheme with boundary conditions at some 
critical high energy scale, such as the grand unification 
(GUT) scale, and keeping the number of free parame- 
ters and assumptions minimal. This is our mSUGRA. In 
this class of models we consider gauge coupling constant 
unification at the GUT scale. In addition, all the mass 
parameters in the supersymmetric breaking sector are 
universal at the same GUT scale. The low-energy sector 
of the model is obtained by evolving all the parameters 
through rcnormalization group equations from the GUT 
scale down to the electroweak scale: this process also in- 
duces the breaking of the electroweak symmetry in a ra- 



diative way. This model is very predictive, since it relies 
only on very few free parameters, but at the same time 
it has a very constrained phenomenology at low-energy. 
It also appears to be quite sensitive to some standard 
model parameters, like the mass of the top and bottom 
quarks (pit and m,b) and the strong coupling constant 
a s . In this class of models there are four free parame- 
ters: the universal gaugino mass parameter Mi/ 2 at the 
GUT scale, the universal soft-scalar mass parameter for 
both the sfermions and the higgscs too at the GUT scale, 
a common trilinear coupling for the third family at the 
GUT scale A and tan (3. The parameters will be varied 
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in the following intervals: 50 GeV < M 1/2 < 1000 GeV, 
< m < 3000 GeV, 1 < tan/3 < 50 and -3 < A < 3. 
The standard model parameters m t , nib and a s are varied 
inside their 2a allowed ranges. 



C. The supersymmetric flux factor T 

The flux factor T defined in Eq. JfjJ) acts as a nor- 
malization factor for the antiproton flux and is a purely 
supersymmetric term. It depends on the mass and cou- 
plings of ncutralinos in the supersymmetric framework 
under study. In Fig. ^ we show the flux factor as a func- 
tion of the neutralino mass for a scan of the eMSSM. 
Fig. |3 reports the case for the mSUGRA scheme. We 
show the values of T separately for the case of comso- 
logically dominant (0.05 < il x h 2 < 0.3, Fig. HJi and 
Figs. 0i, c) and subdominant (Q x h 2 < 0.05, Fig. QJ 3 
and Figs. |2t>,d) relic neutralinos. Among the cosmo- 
logically relevant ones, we also show the configurations 
which yield £l x h 2 inside the preferred range for CDM, as 
determined by the combined WMAP+2dFGRS+Lyman-a 
analysis: 0.095 < n C DAih 2 < 0.131 Q. The results in 
the eMSSM show that the upper values of T are around 
10 -12 GeV -4 for neutralino masses close to the experi- 
mental lower bound (around 50 GeV) and then decrease 
below 10" 14 GeV" 4 for m x ~ 1 TcV. 

In the case of dominant relic neutralinos, the inter- 
val of values for T is restricted, at all masses: in or- 
der to have values of fl x h 2 which fall in the cosmologi- 
cally relevant range, the annihilation cross section inte- 
grated from freeze-out down to the present time must 
be inside the interval 3 • 10 -11 GcV~ 2 < (cr ann w}i nt < 
2 • 10~ 10 GeV~ 2 . We remind that the relic abundance 
depends on (a &nn v) int {VL x h 2 oc (<x a nn«)it)- Tms cross 
section, due to the non-vanishing temperature in the 
early Universe, may differ quite substantially from the 
zero-temperature cross section (a a , nn v)o 1 which is instead 
relevant for the antiproton signal. Usually a correlation 
between (cr ann w}mt and (cr ann w)o is present when the zero- 
temperature (cr ann w)o is large; on the contrary, when 
(fann'i')o is small, temperature corrections in the early 
Universe induce {<7 ann w)mt to deviate, also sizeably, from 
(f ann «}o- This difference in the two cross sections is re- 
sponsible for the band of values of T shown in Fig. [TJi. 

When the neutralino relic abundance is low, such that 
neutralinos are not the dominant component of dark 
matter, T acquires an additional dependence on fl x h 2 
through the rcscaling factor £ 2 . This is shown in Fig. ^d. 
The effect of £ 2 is obviously to reduce T: the lower the 
relic abundance, the smaller £ and then the flux factor. 
The lowest points in Fig. are the ones with lower val- 
ues of fl x h 2 . These configurations, even though they give 
a large (er ann v)o (low fl x h 2 has large values of (cr ann w}i n t, 
and in this case (<J ana v)- m t oc (cr ann z;}o), nevertheless they 
have a low flux factor because they are under-abundant. 
This implies that largely subdominant relic neutralinos 
are likely to provide (almost) undetectable antiproton 



fluxes. This is somewhat at variance with the case of 
direct detection: the difference arises from the fact that 
the antiproton signal (as well as the other Galactic sig- 
nals) depends quadratically on the dark matter density 
(and henceforth on the rescaling factor) while for direct 
detection the dependence is linear and the suppression is 
much milder 0. 

The situation of mSUGRA is shown in Fig. In this 
case, the largest values of the flux factor T are about an 
order of magnitude lower that in the eMSSM: T < 10~ 13 
GeV -4 for light neutralinos. This is a consequence of 
the properties of neutralinos in this constrained type of 
models: neutralinos turn out to be mainly gauginos and 
their couplings, especially to higgs bosons which require 
a mixed higgsino-gaugino neutralino content, are in gen- 
eral smaller than in some sectors of the eMSSM. The 
lower panels in Fig. [5] show the situation in a sector of 
the mSUGRA scheme where the soft scalar masses are 
large: mo > 1 TeV [52, In this sector, the neu- 

tralino may acquire a non-vanishing higgsino component 
|22l I23I ] , as a consequence of the radiative electroweak 
symmetry breaking, and their couplings to higgs bosons 
are enhanced j22|: the consequence on the flux factor is 
in fact a mild enhancement, up to values of T around 
3 - 4 • 10~ 13 GeV -4 , closer to the eMSSM upper values. 



D. The differential antiproton spectrum g(Tp) 

Let us move now to the discussion of differential spec- 
tra of antiprotons which are produced by neutralino an- 
nihilation. The capability of producing antiprotons de- 
pends on the possibility for ncutralinos to produce quarks 
or gluons, cither directly or through decay of their anni- 
hilation products: quarks and gluons will then hadronize 
and eventually produce antiprotons. We have modelled 
the hadronization process by using the PYTHIA Monte 
Carlo [24| ■ The neutralino annihilation is calculated an- 
alytically as described in Ref. j2f|. Neutralino annihila- 
tion occurs at rest in the Galactic frame, and the different 
final states which are open depend therefore on the neu- 
tralino mass. The annihilation may proceed through the 
following channels: production of a fermion pair; produc- 
tion of WW and ZZ; production of a higgs pair; produc- 
tion of a higgs together with a gauge boson (which can 
be the Z boson or the W depending whether the higgs 
is neutral or charged). Apart from the direct production 
of quarks or gluons, the decay chain of the annihilation 
products until a quark is produced is calculated analyt- 
ically. At this stage, the antiproton differential flux is 
obtained from the MC modelling. More details are given 
in App. El 

A sample of p spectra for the four types of neutralino 
annihilation final states is shown in Fig. |3| for different 
values of the neutralino mass: panel (a) shows the spec- 
tra calculated for annihilation into a pure bb state; panel 

(b) refers to annihilation into a pure ZZ state; panel 

(c) refers to an annihilation into a higgs pair, where the 
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FIG. 6: Representative difFerential antiproton spectra per 
annihilation event g(Tp) from neutralino self-annihilation, as 
a function of the antiproton kinetic energy Tp. The different 
curves refer to different neutralino masses: m x = 10 (dotted 
[blue]), 60 (short dashed [black]), 100 (solid [red]), 300 (long 
dashed [green]) and 500 GeV (dot-dashed [magenta]). The 
spectra are selected from our sample of the eMSSM, except 
the one for m x = 10 GeV which refers to an eMSSM with- 
out grand-unification gaugino universality |2?1 WH . All the 
spectra refer to neutralinos with Q, x h 2 = 0.1 and large values 
of the flux factor T. The spectrum for m x = 100 GeV is 
the reference spectrum for the analysis of the astrophysical 
properties in the next Sections. 



scalar higgs has mass of rrih = 120 GeV, the pseudo-scalar 
higgs mass is wia = 200 GeV, for tan/3 = 10 and for a 
vanishing value of the higgs mixing parameter a; panel 
(d) refers to annihilation into a hZ pair, for: m/i = 120 
GeV, tan/3 = 10 and a — 0. Fig. [3]shows the dependence 
of the antiproton spectra on the production energy, fixed 
by the neutralino mass. For instance, in panel (a) the an- 
tiprotons are produced by the hadronization of b quarks 
injected at the energy given by the neutralino mass; in 
panel (b), antiprotons are produced by quarks produced 
by the decay of Z bosons in motion with respect to the 
neutralino rest frame: a Lorentz boost on the hadroniza- 
tion spectra is therefore operative in shifting the fluxes 
to larger kinetic energies. 

Spectra like the ones shown in Fig. [3] are used to cal- 
culate the differential spectra g{Tp). However, as it is 
clear from Eq. (0J, we also need to know the values of 
the branching ratios of each neutralino annihilation fi- 
nal state. The branching ratios will weight the different 
differential spectra, like the ones shown in Fig. [3J An 
example of branching ratios for neutralino annihilation 
is given in Fig. 0] for the eMSSM scheme, and in Fig. 
for the mSUGRA models. 

In the eMSSM, we notice that the annihilation in 
fermions may be sizeable and dominant for masses lower 



than 500 GeV. The two higgses final state is usually of 
the order or lower than the // final state, while the 
gauge bosons final state may dominate, except for very 
large neutralino masses. The mixed gauge+higgs boson 
final state tends to be dominant at very large neutralino 
masses. 

In the case of mSUGRA models, since the neu- 
tralino tends to be a gaugino which couples effectively to 
fermions through sfermion exchange, the // final state 
usually dominates. A relevant production of final states 
other than fermions, especially gauge bosons, occurs in 
the large sfermion masses regime (mo > 1 TcV), where 
sfermion-exchange is suppressed by the large sfermion 
mass and at the same time a higgsino component for the 
neutralino arises: this facilitates both the coupling to 
higgses and to gauge bosons. 

The final result of this analysis is the calculation of 
realistic antiproton differential spectra for neutralino an- 
nihilation. Some representative examples are shown in 
Fig- El f° r different values of the neutralino mass. All the 
spectra refer to neutralinos selected to have Q, x h? = 0.1 
and large values of the flux factor T (close to the up- 
per values of Fig. ^ f° r each mass). All these spectra 
properly take into account all the ingredients discussed 
in this Section: the hadronization spectra and the anni- 
hilation branching ratios. The spectrum for m x = 100 
GeV is the reference spectrum for the analysis of the 
astrophysical properties in the next Sections. The spec- 
tra shown in Fig. are selected from our sample of the 
eMSSM, except the one for m x = 10 GeV which refers to 
an eMSSM without the grand-unification condition be- 
tween the gaugino mass parameters Mi and Mi [2^. |27| . 
In this class of models the neutralino can be as light as 
a few GeV [Mill, m contrast to the standard eMSSM, 
where LEP constraints imply a lower bound on the neu- 
tralino mass of about 50 GeV. For completeness, we have 
therefore included also the representative spectrum for 
m x = 10 GeV, in order to illustrate the effect of propa- 
gation on the primary flux from light neutralinos. How- 
ever, a complete study of the eMSSM without grand- 
unification gaugino universality is beyond the scope of 
this paper. 

Now that we have discussed the source term, we pro- 
ceed to the second step of the calculation: the study 
of how these antiprotons diffuse and propagate in the 
Galaxy and in the Solar System. The result of this analy- 
sis will be the interstellar and top-of-atmosphere (TOA) 
fluxes of primary antiprotons. 



III. DIFFUSION AND PROPAGATION IN THE 
GALAXY 

The propagation of cosmic rays in the Galaxy has 
been considered in the framework of a two-zone diffu- 
sion model, which has been described at length in Rcfs. 
pH IT9L |2&| . Here we only remind the main features of 
this model, and refer to the above-mentioned papers for 
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all the details and motivations. We also present in de- 
tail the quantitative dependence of the secondary and 
primary signal on the propagation parameters. 



The framework 



The disk of the Galaxy is described as a thin disk of 
radius R = 20 kpc, which contains the interstellar gas 
with a surface density £ = 2/misM with h = 100 pc 
and nisM = 1 cm" 3 . It is embedded in a thicker diffu- 
sion halo, supposed to have a cylindrical shape with the 
same radius R as the disk and height L which is not well 
known. The matter density is much lower in the diffusion 
halo so that spallations (rate T = 2/misMC yielding the 
secondary species) of the charged nuclei occur only in the 
disk. Moreover, the standard sources also happen to be 
located in the disk. 

The spatial diffusion of cosmic rays is assumed to occur 
uniformly in the whole (disk and halo) diffusion volume, 
with the same strength. The corresponding diffusion co- 
efficient has been defined as K{E) = K Q P(TZ/1 GV) 5 , 
where 1Z stands for the particle rigidity and Kq and 5 
are free parameters of the model. We also consider the 
possibility that a Galactic wind blows the particles away 
from the disk in the z vertical direction, with a constant 
speed V c . It induces an adiabatic dilution of the energy 
of the particles in the disk due to the sudden change in 
V c . Several other processes modify the antiproton energy 
distribution: ionization losses when interacting with the 
neutral interstellar matter, or from Coulomb losses in a 
completely ionized plasma, dominated by scattering off 
the thermal electrons. To end with, minimal reacceler- 
ation on random hydrodynamic waves, i.e. diffusion in 
momentum space, described by a coefficient Krm related 
to the spatial diffusion K(E), is inevitable |23|. This 
process is assumed to occur only in the disk and is re- 
lated to the velocity of disturbances in the hydrodynam- 
ical plasma Va, called Alfven velocity. In summary, our 
diffusion model has five free parameters Kq, 8, L, V c , Va 
which describe the minimal number of physical effects 
thought to have some role in antiproton propagation. 

The sets of diffusion parameters were constrained in a 
previous work 0] (see also Ref. [3(J) by analysing stable 
nuclei (mainly by fitting the boron to carbon ratio B/C). 
The values we obtained were also shown to be compat- 
ible with the observed secondary antiprotons jlfjj and 
the flux of radioactive isotopes [31|. However, in a first 
step we will disregard these constraints: in order to clar- 
ify which of the propagation parameters are important if 
one wants to compare any possible primary component to 
the background (secondaries), we study the effect of each 
parameter on the signal and the background. Only in a 
second step this additional information on the constraints 
is used to conclude about the variation that will result in 
the primary supcrsymmctric signal (see Sec. lIV^|) . 



B. Solutions for primary and secondary 
antiprotons 



We are interested in the cosmic ray antiproton flux: 
&(r,z,E) = ^-N p (r,z,E) . 

47T 



(8) 



It is related to the differential density N p (r, z, E) = 
dN p (r, z, E)/dE which satisfies the steady-state diffusion 
equation. The general procedure to solve this equation 
as well as references for detailed derivation is given in 
App.[Bj At variance with the solutions already presented 
elsewhere, it proves to be useful, as suggested by the 
study of t he sp atial origin of secondary and primary CRs 
(see Ref. (33 for what is ment by "spatial origin"), to 
introduce the quantities: 



2K{E) 
r - K ^ 

1 sp — 



hT inel (E) - 



(9) 
(10) 



Since many configurations of K(E) — Kq(31Z s and V c 
lead to the same r w and r sp , these new parameters avoid 
automatically a useless discussion about many degener- 
ate values of the diffusion coefficient and make the de- 
pendence over the important parameters more evident in 
formulae (the physical meaning of these new parameters 
is explicited below, see also Ref. [32^. 

The solutions are given below discarding energy re- 
distributions (see App. IB II for the procedure to include 
them). Energetics are not the dominant effects so it is 
interesting to focus on the analytical formulae obtained 
in that case. 

a. The primaries Let us first inspect the primaries: 
the source term is given by <j^ usy (r, z, 2V) described by 
Eq. 0} and discussed in detail in Sec. [n] Primary an- 
tiprotons are produced throughout the whole diffusive 
halo, which is embedded in the dark matter halo. An 
advantage when energy losses and gains are discarded is 
that the solution can be recast as (see App. [BJ) 

N P^ {r = RqiZ = 7 E)= ££Z e (E) x S%Z(R Q , 0, E) 

where the elementary source term (spectrum from a point 
source) given by 

££Z C (E) = Tg(T p ) (11) 
can be separated from the astrophysical part 
00 ( pR 

S*:Z(R Q , 0, E) = ^2 ME, Re) / M&r/R) 
i=i { Jo 

' e -.Av ^ % ~ «)/2] x w (r> x) \ dzrdr . (12) 
sinh^L/2] v >f v ' 

In the above equation, w(r, z) is the effective spatial dis- 
tribution of the primary sources [e.g. p^ M (r,z) for su- 
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persymmetric particles and pdm(^, z) for evaporating pri- 
mordial black holes (PBH)] . We have defined 



ME) 



We also use 



x JoiQRv/R) ■ (13) 



SiL 



A%E) = K(E) 2r-\E) + 2r~\E) + S t coth -I- 



and 



(14) 



(15) 



The functions J and J\ are respectively the Bessel func- 
tions of 0-th and first order, and Q is the i-th zero of Jq. 
The superscript p in A\ (E) indicates that the term r sp 
should be evaluated for the antiproton destruction rate 
and at the parent rigidity. 

Compared to Eq. Q that describes the supersymmet- 
ric source term at each position (r, z) , we isolated in the 
new term £ - P ) "™ e (£') the only required information about 
the production process. For antiprotons produced by 
neutralino annihilations, the flux factor T (see Sec. Ill C|) 
and the elementary spectrum g{Tp) (see Sec. Ill Dl arc 
fully described by the properties of the supcrsymmetric 
and hadronization models. As a result, S , ^ > s r t I ™(i? Q , 0, E) 
is solely dependent on the propagation properties and the 
effective spatial source distribution w(r, z). This function 
is all we need in order to discuss the propagation uncer- 
tainties on the primary fluxes, i.e. the signal detected at 
Solar location Rq. 

b. The secondaries. The secondaries are produced 
from proton sources distributed according to the spa- 
tial supernova remnant distribution in the thin disk 
2hS(z)q(r). These protons are first propagated, leading 
to an equilibrium distribution N p (r, z, E), that in turn 
produces secondary antiprotons when it interacts with 
the interstellar gas. Compared to primaries, secondaries 
diffuse twice. Actually, it is not possible, strictly speak- 
ing, to isolate an elementary source term as for primaries 
(we skip the details, but the interested reader can inspect 
the structure of equations in Ref. 0]). However, it is 
possible to overcome this shortcoming. Antiprotons are 
only produced by protons that are beyond the threshold 
of 7 GeV; in the term A\ that ori gina lly appears in the 
secondary solution (see e.g. Ref. [lj|) ~ and that pre- 
vents from this separation -, one can neglect spallations 
and convection (high energy regime) and approximate 

A\ « K{E) x 2(C«/ii) cothiQL/R) . 

It is then possible to recast the various terms entering 
the solution in order to obtain a formula that (as for pri- 
maries) isolate the dependence on the propagation terms: 



R Q ,0,E) m £iouTce( E ) X B Ztro(R&,Q, E ) 



The corresponding terms are 

Q(E P ) da(E p ,E p ] 



^•source i-Bp) 



dE n 



' Ethtesh K(E p ) dE p 



(17) 



and 



rsoc i r> n pN = V H.i{E,R Q ) 
astro( °' ' ' ~ ^ x 2(C l / J R)coth(C,L/ J R) 

J (Cir/R)2hq disi (r)jrdr . (18) 

This approximate solution is only used to estimate the 
sensitivity of the fluxes to the diffusion parameters. We 
go back to the full solution (see Ref. 01 an d App.lBlfor 
more references) when the final results are presented. 

c. Sensibility to the propagation parameters. With 
the quantities defined above, it is straightforward to eval- 
uate primaries and secondaries fluxes "as if the elemen- 
tary production processes were the same (to focus on 
the astrophysical uncertainties). This defines the rela- 
tive sensitivity to the propagation parameters, and it is 
merely the ratio of the astrophysical part of the signal S 
to the background B: 



Sens\Pax\ = 



(19) 



(16) 



This ratio is likely to depend on the propagation parame- 
ters, in first place because primary sources are located in 
the whole diffusive halo, whereas secondary sources are 
induced spallativcly in the thin disk only. 

We now investigate how the primary flux 5^"™, sec- 
ondary flux -Bistro an d the relative sensitivity Sens de- 
pend on the propagation parameters K(E), r w , r sp , L, 
R, Rq and on the effective source distribution w(r, z). 
This discussion will be general and apply to any primary 
species. It is discussed below for the case of supersym- 
metric primaries, but we will also plot (but not comment) 
the results for PBH antiprotons. 



C. Evolution of fluxes with astrophysical 
parameters 

We now review each one of the above parameters, 
starting with the diffusion coefficient K(E) = Kof31Z s . 
This parameter induces both a change in the normaliza- 
tion - through Kq and only in the high energy regime 
- and in the energy dependence (through 1Z S ). At suffi- 
ciently high energy (above a few tens of GeV), r w , r sp ^ 1 
and Ai and Si become independent on E, so that the sole 
energy dependence 1/K{E) is factored out of Hi(E, Rq), 
i.e. of the Bessel sums. As a result, the quantity 
Sens is insensitive to the choice of K(E), whatever the 
value of the other parameters. There is one subtlety 
left: the secondary elementary production £ ^ lce (E) , 
as defined above, contrarily to the primary's one, is 
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FIG. 7: This plot displays the quantities K{E) x S^jj™ and 
-ff(-E) x -B|st ro , see Eq s - itT^l and JTHt (left panel) and <Sens 
defined by Ea. l|190 (right panel), as a function of the propa- 
gation parameter L (r w = r sp = oo, i.e. no wind, no spal- 
lations) for an isothermal profile. Two cases have been con- 
sidered for the primary signal. The curve labelled supersym- 
metric corresponds to an effective source term proportional 
to w(r,z) = Pumi whereas the curve pbh corresponds to a 
source term proportional to w(r, z) = pdm, as for primordial 
black holes. 



not fully elementary, because it does depend on the 
value of K{Ep) above 7 GeV. However, as we will see 
later, all propagation parameters are designed to have 
about the same K(E) at 100 GeV, so that the quan- 
tity NP> P" m /NP> sec (r = R e , 0, E) is eventually not that 
sensitive to this parameter. 



1. The diffusive halo size L and the radius R of the Galaxy 

These parameters are related to the escape probability 
from the confinement volume (the magnetic halo of the 
Galaxy) . The larger L and R, the greater the probability 
for particles emitted in remote sources to reach us. Ac- 
tually, the side boundary plays almost no role for several 
reasons. First, escape is driven by the closest boundary, 
which is the one at z = ±L as L is likely to be smaller 
than R; second the source distribution is peaked near 
the Galactic center and decreases to very small values at 
large radii (see Ref. [12 for more details). Hence, for 
L < 5 kpc, setting R = 20 kpc or R = oo leaves S and 
B unchanged. The enhancement of fluxes with L can be 
seen in Fig.0(we use here and in other figures the isother- 
mal profile for the dark matter distribution), showing 
Sens[L,r w = oo,r sp = oo] as a function of L (we limit 
the discussion to the supersymmetric case, but the reader 
can straightforwardly conclude for PBH's). For small L, 
only the sources very close to the Solar neighboorhood 
contribute and, as the dark matter source distribution 
is normalized to 1, the supersymmetric and PBH cases 
yield the same value. As L increases, escape is less effi- 
cient and more sources (secondary or primary) effectively 



FIG. 8: Same quantities as in Fig. [7| but as a function of 
the propagation parameter r w i n d, and for three values of L 
(r Bp = oo). 



contribute to the signal. This enhancement is more im- 
portant for primary than for secondaries, as the effective 
number of sources increases respectively as I? (volume 
distribution) and L 2 (surface distribution). In the case 
of primaries, part of the enhancement is also due to the 
mere fact that the number of sources within the diffusive 
box increases with L (we remind that the sources from 
the dark halo to be propagated are those enclosed inside 
the diffusive box, see Ref. Both effects arc respon- 

sible for the evolution of Sens. For L > 5 kpc, no further 
significant enhancement is observed, as the bulk of the 
primary sources (the core radius of the dark matter dis- 
tribution) is then almost entirely enclosed in the diffusive 
halo. To be quantitative, Sens is increased by a factor 
3 for L = 15 kpc compared to L = 1 kpc. Notice that 
the quantity (K(E) x S^™) plotted on the left panel of 
Figs. I7I8I9I does not depend on K(E). To understand 
this property it is sufficient to look at the expression for 



2. The Galactic wind V c throuyh r w 

At high energy (generally a few tens of GeV), propa- 
gation is dominated by diffusion. At low energy convec- 
tion may become the most efficient process (parameter 
r w « 1, see Refs. |32|,|33|) an d it may compete with L 
for escape. The effect of convection is to blow the par- 
ticles away from the disk, leading to an effective size of 
the diffusive halo L* ~ r w . There is a difference with 
the effect of L as the decrease of r w does not lead to 
a decrease of the number of primary sources enclosed 
in the diffusive volume. However, it turns out that the 
effect of the Galactic wind is also more important for pri- 
maries than for secondaries, as the flux is exponentially 
decreased with z for particles created at height z in the 
diffusive halo. This can be seen in Fig. [S] (left panel). We 
clearly see the competition between L and r w in the right 
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FIG. 9: Same quantities as in Fig. [T] but as a function of the 
propagation parameter r sp , and for three values of L (r w = 
oo). 

panel. For large L, the evolution of Sens is completely 
driven by r w , so that we can compare the result to those 
of Fig. [7J When wind is present, the sensitivity to a sig- 
nal is much more reduced than that we would obtained 
with a similar L (i.e. a factor ~ 25 on Sens between 
the case r w = 1 kpc and r w = 15 kpc, compared to a 
factor 3 for L in the same range). For small r w , all pri- 
mary curves converge to the same value, independently 
of L, because then the cosmic rays become blind to this 
boundary, being convected away before having a chance 
to reach the top or bottom of the box. 



3. The relative rate of spallation through r sp 

At low energy, particles can be destroyed more easily, 
because the probability to cross the disk, and thus to 
interact with matter, increases relatively to the escape 
(diffusive or convectivc) probability. The dependence of 
Sens with r sp is displayed in Fig. [21 

When r spa i increases, we are sensitive to sources lo- 
cated farther, and as for L and r w i n( j, the effect is more 
important for primaries than for secondaries. However, 
the effect of r spa i is milder. This is because the cut-off 
due to spallations is less efficient than escape or con- 
vective wind to prevent particles coming from faraway 
sources from reaching us. 



D. A comment about secondaries from GalProp 
model 

Among several other models that are used to describe 
cosmic ray propagation,the fully numerical approach im- 
plemented in GalProp [34| has been widely used. Some 
results obtained within this framework, in particular 
when studying the secondary antiproton spectrum, seem 
to differ (see e.g. Rcf. [35|l from ours, obtained with a 



semi-analytical model. In our paper we want to derive 
constraints on the supersymmetric contribution which 
can be added to the secondary one, when confronting 
with data. Therefore, we take the opportunity of this 
specific work to briefly summarize and discuss some of 
the differences between the two approaches and their re- 
sults. 

First, the approximation that may appear crucial is 
that in order to find analytical expressions for the cos- 
mic ray density, we have to use a simplified description 
of the matter distribution in the Galaxy, whereas with a 
numerical approach, any distribution can be considered. 
However, the results are not strongly affected by this 
hypothesis. In the framework of steady-state diffusion 
models, j3^ has shown that the stable nuclei detected 
in the solar neighborhood were emitted from sources lo- 
cated in a region large enough so that, having sampled 
very different regions of the galactic disc, they are sen- 
sitive to a mean density. Moreover, introducing a radial 
dependence of the matter distribution does not induce 
sizeable difference in the results [3(|. In relation to this 
first point, we have to emphasize that Ref. |35| actually 
does not use a detailed description of the local (i.e. on 
a scale of a few hundreds of pc) gas distribution. As 
a result, they can not provide a reliable analysis of the 
radioactive species, which are very sensitive to the local 
structure of the interstellar medium . 

Second, the numerical approach is still costly in terms 
of computation time, and is less suited to the systematic 
study of different effects. For example, Ref. [34J using a 
predefined small value 6 = 0.3 for the diffusion coefficient 
spectral index, finds that the observed spectrum of B/C 
required small values of the Galactic wind. Indeed, a 
full scan of the parameter space, extended to a range of 
values for 5, revealed that models with higher values of 8 
and with larger values of the galactic wind, were actually 
preferred. This and other results have been thoroughly 
discussed in Ref. j^ij, and also compared with different 
propagation models (such as GalProp). This point is of 
great importance for the present work, as the theoretical 
uncertainties in the antiproton flux are underestimated 
if some parameters are not varied over all their plausible 
values. 

The last relevant difference is actually not related to 
the astrophysical model but to the production cross sec- 
tions. In particular, those relevant for B/C have been 
recently updated (see references in Ref. |35j), whereas 
we use a standard set (see references e.g. in Ref. [l9|)- 
This is a possible way to explain the discrepancy be- 
tween the secondary antiproton flux, but we estimate this 
is unlikely. Indeed, the two sets of cross sections differ 
mainly at low energy, for which the weight of experi- 
mental data is not the greatest. Using the updated set 
should not change the propagation parameters derived 
from B/C and used to propagate antiprotons; the final 
results would essentially remain unaffected. 

To conclude, we do not see any physically relevant dif- 
ference between the two approaches, and they are proba- 
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FIG. 10: The solid lines represent the antiproton flux for a 
m x =100 GeV neutralino and for maximal, median and mini- 
mal astrophysical configurations, for Xb/c — 40. Dotted lines: 
the same, but for for xl/c — 30. The dot-dashed band cor- 
responds to the secondary flux as taken from Ref. O for all 
the configurations giving Xb/c — 40. 

bly equally valuable. There is still some work to be done 
from both sides to understand the origin of the differ- 
ences in the results, which may lie in the methods and 
interpretation of the results, more than in the models 
themselves. 



IV. RESULTS AND UNCERTAINTIES FOR 
THE PRIMARY FLUXES 

We now use all the ingredients previously discussed (as 
well all energy changes) to evaluate the primary inter- 
stellar flux. We try to quantify all the uncertainties that 
could hamper a clear selection or exclusion of supersym- 
metric configurations. They are substantially induced 
by the de gene racy in propagation parameters (see Sec. 

and the choice of a peculiar dark matter 
profile (see Sec. IIV Bl) . 

A. Primary fluxes and related uncertainties 

The propagation parameters have been constrained by 
an analysis of the observed boron over carbon (B/C) ra- 
tio, by means of a Xb/c test over ^6 data points and 
five free parameters [l||. The best Xb/c was f° un d to 
be 25.5. A value of 40 was considered quite conservative, 
corresponding roughly to 4-cr confidence level on B/C 
data interpretation, while a Xb/c = ^0 can be assigned 
to about 2-er confidence level 

In Fig. 1101 we present the result for the primary an- 
tiproton flux for our reference source term for m x = 100 



GeV, whose g(E) is plotted in Fig. We plot the fluxes 
corresponding to the parameters providing the maximal 
and minimal fluxes when all the astrophysical configu- 
rations are taken to be compatible with the analysis on 
stable nuclei, i.e.: Xb/c ^ ^0. For the same set of astro- 
physical parameters we also plot the secondary antipro- 
ton flux. The variation of the astrophysical parameters 
induces a much larger uncertainty on the primary than 
on the secondary flux: in the first case, the uncertainty 
reaches two orders of magnitude for energies Tp < 1 GeV, 
while in the second case it never exceeds 25% (notice that 
these uncertainties are smaller than the nuclear ones, see 
Ref. E3l)- A thorough discussion about why a combi- 
nation of parameters gives the same secondary flux, is 
skipped here but the reader is referred to Ref. [3(j for 
more details. The large variation in the primary signal 
can be understood from the previous discussion: first, the 
exotic signal is more sensitive to astrophysical parameters 
than the standard, as already underlined. Second, this 
has to be weighted by the fact that the secondary flux 
has in its source term an additional K(E). While many 
combinations of K$, <5, L and V c lead to the same sec- 
ondary flux, it is not straigthforward to decipher which 
ones lead to the maximum and minimum primary fluxes. 
Decreasing L and r w decreases the flux, but at the same 
time, to k eep the fit to B/C good, K has to be also 
decreased 0, |3(J , increasing in turn the flux (primaries 
depends on 1/K(E)). However, the first two parameters 
are more important (especially the wind effect) than the 
latter. We give in Tab. [I] the values for these parameters 
yielding the maximum and minimum of the error band in 
both primary and secondary fluxes. The resulting vari- 
ation in Fig. EH can be read off from Figs. I7I8I9I (left 
panels) and Tab. [I] a factor ~ 2000 because of r w and 
i, an additional factor < 4 for r sp (see Fig.0 left panel) 
divided by a factor ~ 50 because of the value of K(E), 
leading to an net scattering of ~ 100. This is almost in- 
dependent on the specific supcrsymmctric configuration. 

As emphasized before, energy redistributions relate 
a specific supersymmetric configuration (by means of 
g{Tp)) to the given propagation configuration. The ef- 
fect on the resulting antiproton flux is expected to be 
mild. We show in Fig. ^2 the result of our analysis for 
the representative eMSSM spectra shown in Fig. cor- 
responding to m x = 60, 100, 300, 500 GeV and for the 
median astrophysical parameters. The low energy be- 
haviour of the fluxes is similar for all the masses: this 
is a consequence of the propagation of the source spec- 
tra, which reduces the intrinsic differences in the origi- 
nal fluxes at low kinetic energies. The high energy be- 
haviour of the fluxes reflects the fact that for higher neu- 
tralino masses, the phase space for antiproton production 
is larger. Since neutralinos in the Galaxy arc highly non- 
relativistic, their mass acts as an effective cut-off on the 
antiproton production kinetic energy 

The effect of propagation on the primary antiproton 
spectrum may also be shown by the following function 
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case 


S 


K 


L 


v c 


V A 


Xb/c 


r w (kpc) 


r sp (kpc) 






(kpc 2 /Myr) 


(kpc) 


(km/sec) 


(km/sec) 


[lGeV/lOGeV] 


[lGeV/lOGeV] 


max 


0.46 


0.0765 


15 


5 


117.6 


39.98 


29./73. 


26./57. 


med 


0.70 


0.0112 


4 


12 


52.9 


25.68 


2.4/9.2 


4.4/15. 


min 


0.85 


0.0016 


1 


13.5 


22.4 


39.02 


0.33/1.8 


0.69/3.1 



TABLE I: Astrophysical parameters giving the maximal, median and minimal supersymmetric antiproton flux and compatible 
wih B/C analysis (xl/c < 40). It is also given in unit of r w , r ap (kpc) for two kinetic energies 1 GeV and 10 GeV. 





FIG. 11: Interstellar primary fluxes calculated as a function 
of the antiproton kinetic energy. The fluxes are calculated 
for the median set of astrophysical parameters. Solid, long 
dashed, short dashed and dotted lines correspond to m x — 60, 
100, 300, 500 GeV, respectively. The fluxes correspond to the 
representative differential antiproton spectra shown in Fig. HU 



FIG. 12: Propagation function CJ£S? of the primary super- 
symmetric antiproton fluxes as a function of the antiproton 
kinetic energy, calculated for the reference fluxes of Figs. 
1611 II Dotted lines refer to m x = 10 GeV, short-dashed to 
m x = 60 GeV, long-dashed to m x = 100 GeV, dot-dashed 
to m x = 300 GeV and solid to m x = 500 GeV. For each set 
of curves, the upper, medium and lower line refer to the max- 
imal, median and minimal set of astrophysical parameters. 
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where &p(0,Tp) is the interstellar antiproton flux after 
propagation, normalized to supersymmetric elementary 
production term. The propagation function C£™v(Tp) is 
a measure of how the source fluxes are deformed by prop- 
agation and diffusion before reaching the solar position 
in the Galaxy and is shown in Fig. ^|for the same rep- 
resentative spectra of Fig. El The energy dependence is 
steeper for low-mass neutralinos, and it becomes some- 
how more symmetric around a maximal values for neu- 
tralinos of increasing mass. The steep rise of C^°^(Tp) 
near the end of the antiproton production phase space 
at Tp = m x is due to reacceleration: while the source 
factor g{Tp) is rapidly vanishing, the propagated flux 
$p(i?Q, 0, Tp) decreases in much milder way because of 
reacceleration effects. This effect is more pronounced for 
the maximal astrophysical configuration, were Va is max- 



imal and it disappears if Va is set to zero. Fig. 1121 also 
shows that the maximal, median and minimal set of as- 
trophysical parameter affect not only the absolute mag- 
nitude of the fluxes but also their energy dependence: 
the distorsion of the original flux differs depending on 
the values of the propagation parameters, as it has been 
discussed in the previous Sections. In particular, the en- 
ergy of maximal transfer for neutralino masses above 60 
GeV shifts from about 1-2 GeV for the maximal set, to 
5-6 GeV for the minimal set. Fig. 1121 shows, at low ki- 
netic energies, a hierarchy in the behaviour of Cg™? (Tp) 
which follows the hierarchy of the neutralino masses: the 
propagation function is larger at low enegies for heavier 
neutralinos, i.e. for harder antiproton fluxes. 

The propagation function C^°^(Tp) can be directly 
used to estimate the propagation effects once the super- 
symmetric production term Tg(Tp) is known. 
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B. Uncertainties related to the dark matter 
distribution 

We have performed all the calculations assuming that 
the Galactic dark matter is distributed as an isothermal 
sphere with a core radius a = 3.5 kpc and local dark 
matter density pi = 0.3 GeV cm~ 3 . For this density pro- 
file, we have estimated that the antiproton propagation 
induces an uncertainty on the primary antiproton flux of 
about two orders of magnitude, especially at low kinetic 
energies. 

Another source of uncertainty on the primary flux 
comes from the shape of the dark matter density pro- 
file, which is only poorly known, and from the allowed 
range of values of pi for any given density distribu- 
tion. We have already commented that for an isothermal 
spherical distribution, the local dark matter density may 
range from 0.18 GeV cm~ 3 to 0.71 GeV cm -3 . More- 
over, the dark matter distribution may be quite different 
from a simple isothermal sphere (see, for instance, Rcfs. 

mmmiH|i3 and referenccs therein): the cold dark 
matter distribution could be non-spherically symmetric, 
it can be singular at the Galactic center, as suggested 
from numerical simulations, or it can present a clumpy 
distribution in addition to a smooth component. Since 
the shape of the Galactic halo enters as Pj iM {r 1 z) in the 
evaluation of the astrophysical part for the primary signal 
•^astro (-^0 ) 0; E), it is a main ingredient in the determina- 
tion of the primary antiproton flux, and the uncertain- 
ties in the description of Pdm(?~, z) may sizeably affect 
the predicted signal. 

The uncertainty in p;, determined by a detailed mod- 
elling of the Galactic component El and mainly due 
to the value of the local rotational velocity [2lJ , depends 
on the shape of the Galactic halo. For the same isother- 
mal sphere, the range in pi may change the primary fluxes 
by a factor which ranges from 0.36 to 5.6: overall, even 
for the simple choice of an isothermal sphere, the antipro- 
ton flux has an uncertainty of a factor of about 15, on top 
of the two orders of magnitude due to antiproton propa- 
gation. We anticipate that, among all the uncertainties 
due to the shape of the Galactic halo, the uncertainty 
coming from pi will turn out to be the most relevant one 
(apart, eventually, the presence of close clumps). 

Independently of the normalization pi, any given den- 
sity profile could, in principle, modify the signal. In 
particular, distribution functions derived from numeri- 
cal simulations are singular towards the Galactic cen- 
ter, where a very high neutralino annihilation rate would 
then occur. We could thus expect that such dark matter 
profiles would induce an enhanced antiproton flux with 
respect to a non-singular distribution. In this class of 
modified density profiles we can also include an isother- 
mal sphere with different values of the core radius a. We 
expect that enlarging the core radius would increase the 
signal. We therefore estimated the modification of the 
cosmic antiproton flux when different core radii and dark 
matter profiles arc used in the source term. The refer- 



ence flux is obtained with our spherical isothermal dis- 
tribution, with core radius a = 3.5 kpc. The results are 
shown in Table |n] Notice that we have used for all the 
profiles the normalization pi = 0.3 GeV cm~ 3 , in order 
to extract the change on the antiproton flux which is due 
entirely to the different shapes of the halos. It is clear 
that each density profile will have to be further imple- 
mented with its own value of pi 0, 13^ . Ecj . From Table 
UHwc notice, first of all, that for small L and r w , we are 
completely blind to what occurs near the Galactic cen- 
ter. Only the very local properties of the dark matter 
distribution are of some relevance for our study. For a 
diffusive halo of 4 kpc, we varied the core radius of the 
isothermal distribution from 2.5 to 5 kpc. With respect 
to our reference values of 3.5 kpc, small a leads to a re- 
duction of the flux by about 20%, while large values of 
a give a 10% increase. For L = 15 kpc - and all the 
other propagation parameters modified consequently - a 
2.5 kpc core radius diminishes the reference flux by 70% 
and a 5 kpc one pushes it up by 25%. The uncertainty 
of a factor two on the core radius of the isothermal dis- 
tribution then reflects in a factor of four indeterminacy 
of the primary antiproton flux. As for the singular den- 
sity profiles, Table |fi] shows that a NFW j3j| distribution 
function does not strongly modify the flux: when com- 
pared to the isothermal case, the flux is increased by no 
more than about 20%, and this occurs when the diffu- 
sive halo size is the largest. For L <5 kpc, the differ- 
ence between an isothermal profile and a NFW singular 
distribution is irrelevant. This result clearly shows that 
it is very improbable for an antiproton produced at the 
Galactic center to reach the Earth. 



Finally, one can deal with halos which contain regions 
of enhanced density called clumps. In this sub-halos, the 
neutralino annihilation is more effective and the signal 
can be increased by some enhancement factor. However, 
as also suggested by Ref. ^lj, this enhancement is not 
propagation dependent and simply acts on the antiproton 
flux as a normalization factor. From the analysis of Ref. 
|4l| , the average enhancement is likely to be smaller than 
a factor of 5. A detailed analysis of this point is beyond 
the scope of our paper; however, the effects of such an 
enhancement are briefly discussed at the end of Sec. El 



In conclusion, we wish to remark that our choice of an 
isothermal sphere with a core radius a = 3.5 kpc and local 
dark matter density pi = 0.3 GeV cm -3 , together with 
the best choice for the astrophysical parameters which 
govern diffusion and propagation in the Galaxy, represent 
an optimal choice for the prediction of the antiproton 
signal. Our results will not be dramatically modified by 
a different choice for the density profile, while a different 
choice for the local dark matter density is easily taken 
into account as a normalization factor. 
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L (kpc), TV, r sp 


Qprtin Qprim 
^a = 2.5 _ ^ref 


Qprim Qprtin 
^ a =5 ~ D re/ 


NFW ~ °ref 


Qprim 
a re.f 




Qprim 


15, 28.66, 25.54 


-69.5% 


+23.9% 


+ 19% 


4, 2.38, 4.41 


-21.5% 


+9.9% 


-0% 


1, 0.33, 0.69 


< 1% 


< 0.2% 


-0% 



TABLE II: Sensitivity to the core radius of an isothermal profile, and comparison of the NFW and istothermal profiles, for 
three representative propagation sets at Tp — 1 GeV. These propagation parameters correspond to the minimum, median and 
maximum primary flux compatible with nuclei analysis. The reference value S^™™ is for an isothermal halo whose core radius 
is a — 3.5 kpc. Notice that for higher energies, the results would be the same than those provided by the set L = 15 kpc (purely 
diffusive transport). 




FIG. 13: Top-of-atmosphere antiproton fluxes as a function 
of the antiproton kinetic energy for the m x = 100 GeV refer- 
ence case. The upper (lower) set of curves refer to the maxi- 
mal (minimal) set of of astrophysical parameters. Solid curves 
show the interstellar fluxes. Broken curves show the effect of 
solar modulation at different periods of solar activity: (j) = 500 
MV (long dashed), <f> = 700 MV (short dashed), <f> = 1300 MV 
(dotted). 



V. TOP-OF-ATMOSPHERE FLUXES: 
COMPARISON WITH DATA AND RESULTS 
FOR SUPERSYMMETRIC MODELS 



Now that we have calculated the interstellar fluxes 
of antiprotons at the Sun's position in the Galaxy, we 
have to further propagate them inside the heliospherc, 
where the cosmic ray particles which eventually reach 
the Earth are affected by the presence of the solar wind. 
We model the effect of solar modulation by adopting the 
force field approximation of the full transport equation 
[42^ . In this model, the top-of-atmosphere (TOA) an- 
tiproton flux $J OA is obtained as: 



<j>J OA (ffJ OA ) 



„TOA 



,,1S 



FIG. 14: Primary TOA antiproton fluxes as a function of the 
antiproton kinetic energy, for the representative spectra of 
Figs, ©in the eMSSM. The solid line refers to m x = 60 GeV, 
the long-dashed line to m x = 100 GeV, the short-dashed 
line to m x = 300 GeV and the dotted line to m x = 500 
GeV. The astrophysical parameters correspond to the median 
choice. Solar modulation is calculated for a period of minimal 
solar activity. The upper dot-dashed curve corresponds to 
the antiproton secondary flux taken from Refs. Ugl. l28l . Full 
circles show the BESS 1995-97 data |4*H| : the open squares show 
the BESS 1998 data 0| ; the stars show the ams data and 
the empty circles show the caprice data 



where E and p denote the total energies and momenta of 
interstellar and TOA antiprotons, which are related by 
the energy shift: 



El s 



(22) 



(21) 



where the parameter <f> is determined by fits on cosmic 
ray data. In our analysis, we will adopt the value <p = 500 
MV for periods of minimal solar activity, corresponding 
to the years around 1995-1998, = 700 and <j> = 1300 
MV for a transient period and for the solar maximum, 
respectively, which will be used for years 1999 and 2000. 

Fig. E| shows the TOA antiproton fluxes for the 
m x = 100 GeV reference configuration and for the max- 
imal and minimal sets of astrophysical parameters. The 
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Ti° A (GeV) 

FIG. 15: Primary TOA antiproton fluxes at solar maximum 
for the transient periods of solar activity of years 1999 and 
2000. The upper set of curves show the antiproton secondary 
fluxes. The lower set of curves show the primary antiproton 
fluxes obtained for the representative m x = 100 GeV case. 
The solar modulation parameter is fixed at 700 MV (solid 
lines) and at 1300 MV (dotted lines). The astrophysical pa- 
rameters correspond to the median case. Stars and full circles 
correspond to BESS 1999 and 2000, respectively EBTl . 



figure shows that solar modulation has the effect of de- 
pleting the low-energy tail of the antiproton flux. The 
effect is clearly more pronounced for periods of stong so- 
lar activity, when the solar wind is stronger. 

Data on antiprotons at Earth are now abundant, 
mostly after the missions of the balloon borne detector 
BESS. This experiment has provided measurements at 
different periods of the solar activity [43l \aA |4o| . It has 
now collected more than two thousand antiprotons be- 
tween 200 McV and 4 GeV. Data at solar minimum have 
been taken also by the AMS experiment on board of Shut- 
tle 0| in an energy range similar to BESS, and by the 
caprice balloon at higher energies, namely between 5 
GeV and 40 GeV 0- All the data at solar minimum 
are plotted in Fig. ^] along with the secondary refer- 
ence flux (for details, see Ref. 0) and our predictions 
for primary fluxes at different neutralino masses in the 
eMSSM: m x = 60, 100, 300, 500 GeV and for the median 
set of astrophysical parameters. We notice that the pri- 
mary flux from neutralino annihilation is at most of the 
same order of magnitude as the secondary flux, and this 
occurs for neutralino masses close to their current lower 
bound in the eMSSM, which is around m x ~ 50 GeV. 
We remind that the representative supcrsymmctric con- 
figurations plotted in Fig. 1141 refer to a large antiproton 
production for each mass (i.e. they correspond to large 
values of the T parameter shown in Fig. This indi- 
cates that the antiproton signal for neutralino dark mat- 
ter will hardly produce an excess over the secondary flux, 
for an isothermal matter profile of the Galactic halo and 
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FIG. 16: Antiproton flux at solar minimum from neutralino 
annihilation calculated at Tp = 0.23 GeV, as a function of 
the neutralino mass for a generic scan of the eMSSM. The 
flux is calculated for a smooth halo described by an isother- 
mal profile with core radius a = 3.5 kpc and for the me- 
dian set of astrophysical parameters. Crosses (in red) refer 
to cosmologically dominant neutralinos (0.05 <Q, x h 2 < 0.3); 
dots (in blue) refer to subdominant relic neutralinos (Q x h 2 < 
0.05); light circles (in green) show the eMSSM configura- 
tions for which the neutralino relic abundance lies in the 
preferred range for CDM, as determined by the combined 
WMAP+2dFGRS+Lyman-a analysis: 0.095 < SlcDAih 2 < 
0.131 Q. The shaded region (in yellow) denotes the amount 
of antiprotons, in excess of the secondary component [lH [. 
which can be accomodated at Tp = 0.23 GeV in order not to 
exceed the observed flux, as measured by BESS (43l. l44j . All 
the points of the scatter plot which lie below the horizontal 
black line are compatible with observations. 



for pi = 0.3 GeV cm -3 . This occurs for the median (and 
best) choice of the astrophysical parameters which gov- 
ern the diffusion and propagation of antiprotons in the 
Galaxy. Clearly, the maximal set of astrophysical pa- 
rameters, which produces fluxes about one order of mag- 
nitude larger than the median set, may produce a large 
excess, for neutralino masses below 100-200 GeV. This 
excess could then be used to constrain supersymmetric 
models since the secondary flux is perfectly compatible 
with the data. However, for setting constraints on su- 
persymmetry in a conservative way, we should instead 
use the set of astrophysical parameters which produces 
the minimal fluxes. In this case, the primary fluxes are 
lower than the ones plotted in Fig. E]by about one or- 
der of magnitude, as discussed in the previous Section. 
In conclusion, our analysis shows that, due to the large 
uncertainties in the primary fluxes, the antiproton signal 
is not suitable at present for setting conservative con- 
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FIG. 17: Antiproton flux at solar minimum from neutralino 
annihilation calculated at Tp = 37.5 GeV, as a function of the 
neutralino mass for a generic scan of the eMSSM. Notations 
are as m Fig. ED The shaded region (in yellow) denotes the 
amount of antiprotons which would be required at Tp — 37.5 
GeV in order to explain the possible excess in the BESS data 
|43ll44l over the secondary component Hsl . All the points of 
the scatter plot which lie below the upper horizonthal black 
line are compatible with observations. 



straints on supcrsymmctric models. For this wc need a 
better knowledge of the astrophysical parameters that 
govern the diffusion and propagation of primary antipro- 
tons in the Galaxy. 

Antiproton data are also available for periods of intense 
solar activity from the BESS detector. Fig. H5lshows these 
data together with the secondary flux and the primary 
flux calculated for the representative m x = 100 GeV 
configuration. The astrophysical parameters are fixed 
at their median values. Also at solar maximum, we see 
that the secondary flux is compatible with the data and 
the supersymmetric flux is sensitively smaller than the 
secondary one. 

In our discussion so far, wc have commented that con- 
servative and solid constraints on supersymmetric models 
require the use of the minimal set of astrophysical param- 
eters. This attitude is needed in setting limits. However, 
the best and most probable choice of astrophysical pa- 
rameters is the median one, and we will therefore adopt 
from here on this set of parameters for our analyses. It 
is likely that a sharpening of the knowledge of the propa- 
gation parameters will lead to a shrinking of the allowed 
uncertainty band around the central (median) value. 

In order to compare the experimental results with a 
full scan of the supersymmetric parameter space, we cal- 
culate the TOA antiproton fluxes in two different energy 
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FIG. 18: The same as in Fig. 1161 for a generic scan of the 
mSUGRA scheme. 



bins and compare our results with the excess which can 
be accomodated above the secondary flux in order not 
to enter in conflict with the experimental data in that 
energy bin. We have chosen a low energy bin: Tp = 0.23 
GeV, and a high energy one: Tp = 37.5 GeV. As can be 
seen in Fig. 1141 in the low energy bin the secondary flux is 
perfectly compatible with the data, therefore no excess is 
needed: this allows us to set an upper bound on the possi- 
ble amount of antiprotons of primary origin which can be 
accomodated: $J OA (Tp = 0.23 GeV) < 2.09 • 10 -3 m~ 2 
s _1 sr" 1 GeV -1 . This value is obtained by taking into 
account the values and uncertainties of both data and 
secondary flux at T p = 0.23 GeV. At T p = 37.5 GeV, 
even though the data and the secondary flux are statisti- 
cally compatible, a possible excess may be accomodated, 
since the central value of the experimental point indicates 
a much larger flux as compared to the secondary compo- 
nent. In this case, we can define an interval of values for 
a possible excess: 0.04 ■ 10~ 3 < $J OA (Tp = 37.5 GeV) < 
1.87- 10 -3 in units of: m -2 s _1 sr" 1 GeV -1 . We compare 
these intervals with our calculations in the eMSSM and 
in mSUGRA. 

Fig. shows the scatter plot of the antiproton flux 
calculated at Tp = 0.23 GeV for a generic scan of the 
eMSSM scheme. The supersymmetric fluxes are clearly 
largest at low neutralino masses, and they fall down as 
the neutralino mass increases mostly because the neu- 



tralino number density in the Galaxy scales as m~ . A 
small fraction of configurations with masses below 100 
GeV can provide fluxes which could be potentially at the 
edge of producing an excess, but wc remind here that for 
a safe exclusion of these configurations wc should use the 
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FIG. 19: Correlation between the antiproton flux at Tp = 0.23 GeV and Tp = 37.5 GeV shown in Figs. 1161 and Il7l for the 
astrophysical enhancement parameter r\ = 3 and 10 (overdense halos) and for the median set of astrophysical parameters. 
Circles (in green), dots (in blue) and crosses (in red) denote configurations with neutralino masses in the ranges: 50 GeV 
< m x < 150 GeV, 150 GeV < m x < 300 GeV and 300 GeV < m x < 1 TeV, respectively. The horizonthal line denotes 
the upper limit on the antiproton flux at Tp = 0.23 GeV coming from the BESS data [igl. |44| . once the secondary component 
[13 is taken into account. The rightmost vertical line denotes the corresponding upper limit at Tp = 37.5 GeV. The shaded 
area indicates configurations which can explain the possible excess in the data |43l. |44| over the secondary component [Tsj at 
Tp = 37.5 GeV, without giving an excees at low kinetic energies (Tp — 0.23 GeV). 



minimal set of astrophysical parameters, which provides 
a flux which is about one order of magnitude smaller. In 
any case, a reduction of the uncertainties on the primary 
flux calculation and a future reduction of experimental 
errors may eventually either allow to set limits to super- 
symmetry or show a positive excess of antiprotons in this 
low-energy bin, a fact which could then be explained as 
originated by ncutralinos of masses below 100 GeV. Fig. 
ITTl shows the scatter plot of the antiproton flux calcu- 
lated at T p = 37.5 GeV for the same scan of the eMSSM. 
In this case we observe that all the supersymmetric con- 
figurations are compatible with data, but there are no 
supersymmetric models which can allow us to explain 
the discrepancy between the data and the secondary flux 
as due to an excess of supersymmetric origin. 

The situation in the mSUGRA scheme is shown in Fig. 
^] with the flux of antiprotons at Tp = 0.23 GeV. In this 
case, as already observed in connection with the proper- 
ties of the mSUGRA source term, the antiproton fluxes 
are smaller than in the eMSSM case. Nevertheless, a re- 
stricted fraction of mSUGRA configuration is potentially 
cxplorablc in the future, with a reduction of the experi- 
mental error of about a factor of 2-3. 

The fluxes we have shown so far, all refer to a dark mat- 



ter density distribution in the form of an cored isothermal 
sphere. Clearly a halo profile which is able to produce an 
overdensity with respect to the isothermal sphere would 
produce a larger antiproton flux. We can parameterize 
the enhancement of the matter density by a multiplica- 
tive factor rj, which then enters as rj 1 in the calculation 
of the antiproton primary flux, since the flux depends 
on the square of the matter density. The origin of the 
overdensity may be due, for instance, to flattening of the 
Galactic halo or to the presence of clumps. It the latter 
case, the enhancement factor is likely to be smaller than 
about 5, once the result of Ref. are implemented with 
our discussion on the antiproton diffusion region in the 
Galaxy. The enhancement factor may also be related to 
a different choice of the local dark matter density, which 
has been fixed at pi =0.3 GeV cm ~ 3 in our analysis. We 
remind that our cored isothermal sphere allows factors of 
77 up to: rj - (0.71/0.3) = 2.4 [2l|. Clearly, a complete re- 
analysis of the propagation and diffusion properties will 
be required for each different choice of the halo shape: 
this reanalysis will give the amount of enhancement con- 
cerning the specific halo. Anyhow, regardless of how the 
enhancement 77 is obtained, we can discuss the effect of 
such an increased flux by using rj as a normalization fac- 
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tor to show the amount of enhancement which is required 
in order to interpret the antiproton excess at Tp = 37.5 
GeV as due to neutralino dark matter, without exceeding 
the upper limit on the antiproton flux at Tp = 0.23 GeV. 
Fig. ^)] shows the correlation between the eMSSM an- 
tiproton fluxes at T p = 0.23 GeV and Tp = 37.5 GeV for 
77=3 and 10. The supcrsymmetric configurations which 
could fulfill this requirement are the ones which fall in- 
side the shaded area in Fig. ^| We see that the possible 
excess at T p = 37.5 GeV requires halo overdensities of at 
least a factor of 2-3 and neutralino masses larger than a 
few hundreds of GeV. This last property is simply un- 
derstood on the basis of the fluxes shown in Fig. [21 
the phase space cut off at Tp = m x implies that light 
neutralino would need a huge overdensity factor in or- 
der to match the observed antiproton flux at Tp = 37.5 
GeV, but this would produce an exceedingly large flux at 
Tp = 0.23 GeV. On the contrary, heavy neutralinos have 
a phase space cut off at much higher kinetic energies, 
and therefore a mild overdensity may enhance the flux 
at Tp = 37.5 GeV without giving conflict at low kinetic 
energies. 



VI. CONCLUSIONS AND PERSPECTIVES 

We have calculated the flux of antiprotons produced 
by relic neutralino annihilations in the Galactic halo in 
a detailed diffusion model constrained by analysis of sta- 
ble and radioactive nuclei. The source of antiprotons is 
studied both in a low-energy minimal supcrsymmetric 
standard model (eMSSM) and in a supcrgravity inspired 
supersymmetric scheme (mSUGRA). We find that the 
interstellar primary antiproton fluxes are affected by a 
large uncertainty, which spans two orders of magnitude 
at low antiproton kinetic energies. This is at variance 
with the secondary antiproton fluxes (whose uncertainty 
never exceeds 24% O) and it is mainly related to the 
fact that the source of the primary flux is located inside 
the diffusive halo, whose size is unknown. By adopting 
a conservative choice for the dark matter density-profile 
and propagation parameters, no supersymmetric config- 
uration can be excluded on the basis of an excess over 
the existing data. Actually, the data are quite well ex- 
plained by the secondary contribution alone. However, 
if we adopt the best values for the propagation param- 
eters (corresponding to a thickness of the diffusive halo 
of 4 kpc), a window of low-mass neutralino configura- 
tions provides fluxes which, once summed up to the sec- 
ondary contribution, are in excess of the experimental 
measurements. We have shown that the sensitivity to 
the antiproton signal is increased with the halo size and 
limited by strong convection. An improved knowledge of 
the propagation parameters will certainly help in reduc- 
ing the uncertainty on the primary flux and consequently 
it could allow us to set more severe constraints on the su- 
persymmetric parameter space. 

The sensititity of the primary antiproton flux on the 



shape of the dark matter density profile has also been 
investigated. We have shown that the shape of the dark 
matter density distribution docs not introduce large un- 
certainties. In particular, we have demonstrated that a 
NFW distribution can increase the primary antiproton 
flux by no more than 20% with respect to the isother- 
mal profile. Indeed, it is very improbable to detect at 
Earth antiprotons produced in the centeral regions of the 
Galaxy, where the two distributions differ most. 

In the next years several balloon-borne experiments 
such as BESS, space-based detectors such as AMS and 
satellites as PAMELA will provide very abundant and ac- 
curate data on the antiproton flux. At the same time, 
new data on cosmic ray nuclei are expected and would 
lead to a better knowledge of the cosmic-ray propagation 
mechanisms. We could thus expect a dramatic reduction 
of the uncertainties affecting the neutralino-induced an- 
tiproton flux and much more definite predictions for an- 
tiprotons of supersymmetric origin will then be possible. 
Many efforts are also devoted to other indirect neutralino 
searches, such as positrons and antidcutcrons in cosmic 
rays, gamma rays and up-going muons, as well as to di- 
rect searches in deep underground laboratories, giving 
thus the hope that more constraining analysis on the ex- 
istence of relic neutralinos in the halo of our Galaxy will 
be viable. 



APPENDIX A: CALCULATION OF THE 
ANTIPROTON DIFFERENTIAL SPECTRUM 
PER ANNIHILATION EVENT 

The antiproton differential spectrum per annihilation 
event g{Tp) is calculated by following analytically the de- 
cay chain of the neutralino annihilation products until a 
quark or a gluon h is produced. The antiproton spec- 
trum is then obtained by a Monte Carlo modelling of 
the quark and gluon hadronization (we make use of the 
PYTHIA package (24|). We have produced the p differ- 
ential distributions for h = u, d, s, c, 6, gluon at various 
injection energies for each h (the t quark is assumed to 
decay before hadronization and is treated analytically, 
since in addition to its standard model decay into W + b, 
it may have a supersymmetric decay into H + b). When- 
ever we need the p distribution for an injection energy 
different from the produced ones, we perform an inter- 
polation. In order to obtain the antiproton differential 
distribution in the neutralino rest frame we perform the 
necessary boosts on the MC spectra. 

For instance, let us consider a p production from a 
neutralino decay chain of this type: 

XX -—>A—>a^h-^p . (Al) 

The antiproton differential spectrum per annihilation 
event g{Tp) is then obtained by the product of the 
branching ratios for the production of A, a and h in the 
decay chain, with the differential distribution of antipro- 
tons produced by the hadronization of an h injected at an 
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energy -E pro d (denned in the rest frame of the a decaying 
particle) double boosted to the x reference frame: 



g(T p ) = BR(xx A)BR(A -> a)BR(a h) x 



dTf, 



boost a — >A. 



(A2) 



boost A — *x 



The first boost transforms the spectrum from the rest 
frame of a (in which h is injected with energy -E pro d) to 
the rest frame of A. The second brings the distribution 
to the rest frame of %■ Each boost is obtained by the 
following expression: 



g(Ep) = 77 



E' 



dN% 
~dE ! 



dE' 

7/3 p' 



(A3) 



where Ep = Tp + mp is the total antiproton energy, p' = 
(E^ + rrip), 7 and f3 are the Lorentz factors of the 
boost and the interval of integration is defined by: 



E' ± 



E pmd ,"fEp 1±/?a/1 



(A4) 



APPENDIX B: SOLUTION OF THE DIFFUSION 
EQUATION 

In cylindrical geometry, the differential density 
N p (r, z,E) is given by 



= 



K{E) 



d 2 







-Vc- 



d 



dz 2 r dr \ dr J J ^ dz 
- Qf{r lZl E)-2h5{z)T p {E)N p {r lZ ,E) , 



N p {r,z,E) 
(Bl) 



where the energy losses have been omitted, for the sake 
of clarity. The source term includes primary antiprotons 
- from exotic sources present in the dark halo, annihi- 
lating throughout the diffusive halo of the Galaxy -, and 
secondary antiprotons - standard p and He CRs spallat- 
ing on the interstellar gas in the thin disk -, and may be 
written as 

Q p (r, z, E) = q p - plim (r, z, E) + 2hS(z)q p - scc (r, 0, E) . 

(B2) 

A convenient way to solve Eq. (|B1|I is to expand all 
the functions f(r) (the density N(r) and the source dis- 
tribution q(r)) that depend on r on the orthogonal set of 



Bessel functions { Jo(Qx)Y~ (where Q are the zeros 
of Jq). These Bessel transforms arc defined as 



(B3) 



with 



fi 



j2ir\ I P-f(P R ) J o(CiP)dp 
J i IW J o 
Using the Fourier-Besscl 



(B4) 

N p {z,E) and 



coefficients 

qf(z,E), there is no conceptual difficulty to extract so- 
lutions of Eq. (|Bip . We do not repeat the derivation 
that can be a bit lenghty. Solutions for primaries can be 
found in Ref. 01, whereas the one for secondaries has 
been discussed in Ref. . 



1. Energy losses, tertiary component 

Following the procedure described e.g. in Ref. [T^ |. 
energy losses and diffusive reacceleration lead to a differ- 
ential equation on Ni(z = 0, E) 



AiNi(0,E) = Qi(E) 



(B5) 



where bi, 



+ b 



Coul 



badiab includes the three 



kinds of energy losses. The exact forms for these terms 
may be found in Refs. 0, |2^. The resolution of this 
equation proceeds as described in App. (A. 2), (A. 3) and 
App. (B.) of Ref. to which we refer for further 

details. The source term also takes into account the so- 
called tertiary component, q\ er {E), corresponding to in- 
clastic but non-annihilating reaction of p on interstellar 
matter. This mechanism merely redistributes antipro- 
tons towards lower energies and tends to flatten their 
spectrum. 
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